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Abstract 

We establish a probabilistic representation for a wide class of linear deterministic p.d.e.s with 
potential term, including the wave equation in spatial dimensions 1 to 3. Our representation 
applies to the heat equation, where it is related to the classical Feynman-Kac formula, as well 
as to the telegraph and beam equations. If the potential is a (random) spatially homogeneous 
Gaussian noise, then this formula leads to an expression for the moments of the solution. 
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1 Introduction 



The purpose of this paper is to present a form of the Feynman-Kac formula which apphes to a wide 
class of hnear partial differential equations with a potential term, and, in particular, to the wave 
equation in dimensions d < 3. In the case of the heat equation, this gives an expression that differs 
from the classical Feynman-Kac formula. As an application, we consider a random potential term 
which is a spatially homogeneous Gaussian random field that is white in time. In this case, our 
approach provides a probabilistic representation for all product moments of the solution, which has 
already shown its usefulness (see [TO]). 

We begin by giving an informal derivation of the representation in the special case of the heat 
equation with potential, where we can contrast it with the classical Feynman-Kac formula. Consider 
the heat equation on M*^ with a deterministic potential V{t,x): 



du(t, x) 

dt 
u{0, x) 



-Au{t, x) + V{t, x)u{t, x) 
fix). 



:i.v. 



The classical Feynman-Kac formula for the solution {u{t,x), t > 0, x G 
states that, under appropriate conditions on V and /, 



(see for instance 



u(t, x) 



E. 



B 



/(5i)exp^ V{t-s,B,)ds 

and is the expectation for Brownian motion 



where {Bt, t > 0) is a Brownian motion in 
started at Bq = x. 

We now heuristically derive an alternative probabilistic representation to (11. ip . which will be 
rigorously established as a special case of the main result in section [31 We start by writing Duhamel's 
formula for the solution u{t,x), using the Green's function, as follows: 



u(t, x) 



Pt{x - y)f{y)dy + 



JM.'^ 



Pt-s{x - y)V{s, y)u{s, y)dyds, 



'1.21 



where 



Pt[x) 



exp 



\x\ 



2t 



(27rt)'^/2 

We use (11.21) as the start of an iteration scheme. Substituting this expression for u{s,y) back into 
the right hand side of ( II. 2p suggests the following series expansion for u{t, x): 



u 



{t,x) = ^Im{t,x), 



1.3) 



m=0 



where Jo(t, x) = J^^Ptix - y)f{y)dy and 

Im+lit, x) 



Jw^ 



Pt-s{x - y)V{s, y)Im{s, y)dyds. 



;i-4) 



We wish to write an explicit expression for Im{t,x). To begin with, let 

w{t, x) = /o(t, x). 
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For convenience, let Sn+i = t and i/n+i = x. Then we have 
where 



/ / \T\Psk+i-Skiyk+i-yk)V{sk,yk)]w{si,yi)dyds 



;i.5) 



Tmit) = {(si, . . . , Sm) : < Si < ■ ■ ■ < < i:}, 

dy = dyi ■ ■ -dym and ds = dsi ■ ■ -dsm- An alternative derivation of the series expansion fll.3p for 
u{t,x) starts by expanding the exponential in the classical Feynman-Kac formula as a Taylor series 
and confirming that the terms correspond to the expansion (11 .Sp . However, we will make use of the 
iterative formula (ll.4p later on. 

A basic observation is that domain of integration Tm{t) has volume t™/m!, which, except for 
a missing exponential factor, is a Poisson probability. If N{t) is a rate one Poisson process, then 
P[N[t) = m] = t^e'* /m\. Let Ti < T2 < ■ ■ ■ be the times of the successive jumps of the Poisson 
process, and let tq = 0. It is well known that if we condition on Nf = m, then the vector (ri, . . . , r^) 
is uniformly distributed over the simplex Tmit). The time reversed sequence t — r^, . . . , t — ti is also 
uniformly distributed on Tm(t). Therefore, setting Sk = t — r,n+i~k and replacing y^ by y^+i-k, so 
that yo = X, we may rewrite the expression (11.51) for Imit, x) as 



/ YlPrk-n-iiVk - yk-l)V{t - Tk,yk) W{t - T^,ym)dy l{Nit)=m} 



where is the expectation with respect to the Poisson process. But we can also exploit the fact 
that pt{x) is the probability density for the increments of a c?- dimensional Brownian motion. Thus 



/ YlPrk~Tk-i{yk - yk~i)V{t - Tk,yk) w{t - Tm,ym)dy 



E, 



B 



\[V{t-Tk,B^^)\w{t-T„,,BrJ 



\k=l 



where E^ denotes the expectation with respect to Brownian motion started at x, and therefore. 



Summing over m, we get 

M(t, x) = e'E^E^ 



' Nt \ 



\k=l 



;i.6) 



The representation (11.60 . unlike the classical Feynman-Kac formula, does not use the entire 
Brownian path but only the values at a finite (random) set of times. This allows us to extend this 
type of representation to equations where the differential operator is not the infinitesimal generator 
of a Markov process. All we will require is a Poisson process and an independent stochastic process 
whose one dimensional marginals give the Green's function for the differential operator. In particular, 
we will treat the case of the wave equation with potential in dimensions d < 3. 
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The outline of this paper is as follows. In Section [2], we describe the class of equations that we 
will consider and establish a series representation in the case of a bounded potential. In Section [31 
we establish our Feynman-Kac-type formula analogous to fll.6p . where the Brownian motion will be 
replaced by a suitable spatial motion that depends on the particular equation being considered. In 
Section |U we give an application to the situation where the potential is a Gaussian random field 
whose covariance is formally given by 



E 



F{t,x)Fis,y) =Soit-s)fix-y). 



In this equation, S{-) denotes the Dirac delta function, / : M'^ ^ M is continuous on M'^ \ {0} and 
the right-hand side is such that f{x — y) is indeed a covariance. This type of covariance is widely 
used in the literature, including for instance O [6l [H [T7j. In this case, unlike the classical Feynman- 
Kac formula, the noise is too rough for the probabilistic representation of the solution to make 
sense. Instead, we establish in this section a formula for the second moment of the solution. Section 
[5] contains the extension to ra-fold product moments. This formula makes use of a Poisson random 
measure combined with a spatial motion. The first two named authors have made use of this formula 
[To] to establish intermittency properties of the solution to the wave equation with potential. 

We end by making a few comments on related literature. Probabilistic representations of the 
solution to deterministic p.d.e.'s abound. The closest related work seems to be results on random 
evolutions, surveyed for example in Hersch [12] and Pinsky [TH]. These references give probabilistic 
representations for some hyperbolic equations, including the Poisson representation for the damped 
wave equation in one spatial dimension (also known as the telegraph equation) as developed by 
Marc Kac [T3l [T4] . Related also is the use of random flight models for the Boltzmann equation, as 
described, for example, in [191 HI]. We cannot quite find our approach represented in this literature. 
The use of Poisson probabilities is also implicit in other works: Albeverio and Hoegh-Krohn [Ij 
and Albeverio, Blanchard, Coombe, Hoegh-Krohn, and Sirugue [2] have used the idea that the 
multiple integrals involved in the expansion of Feynman integrals are related to Poisson probabilities 
(see also [23]). AH these works display the usefulness of probabilistic representations in studying 
problems of asymptotics, homogenization and perturbation theory for the deterministic p.d.e. For 
parabolic equations with random potentials, the classical Feynman-Kac formula has been a key tool, 
for example in the parabolic Anderson problem (see Carmona and Molchanov [3]) and in random 
waves (see Oksendal, Vage and Zhao [H]). We hope that our representation may play a similar role 
for other equations with a random potential. 

2 Series representation for bounded potential 

Our probabilistic representation will be for the integral equation 

u{t, x) = w{t, x) + ds S{s, dy)V{t — s, X — y)u{t — s, X — y). (2.1) 
Jo Jr'' 

In this section, the key assumption is the following: 

Assumption A. For each t > 0, S{t, dy) is a signed measure on Mf^ satisfying 

sup |5(t,M'^)| < oo for all T > 0, 

te[o,T] 

where |S'(t,M'^)| denotes the total variation. 
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It is well known that a large class of linear partial differential equations of the form 

Lu(t,x) = V{t,x)u(t,x) 

can be recast, using their Green's functions, into this integral form. We briefly recall some illustrative 
examples that we consider later. 

Example 2.1. (a) The heat equation on W''. Take L = ^ — |A and S{t,dy) = pt{y)dy. Then 
for a suitable initial condition u{0,x) = fo{x), the Green's function representation of the heat 
equation leads to the integral equation (12. ip with 



w{t, x) 



foix - y)pt{y)dy. 



(b) The wave equation on for d < 3. Take L = ^ — A and 



r 1 



S{t,dy) 



^^{\y\<t}dy if (i = 1, 

^ l{lyl^t}dy ifd = 2, 

ff d = 3, 



4TTt 



where crp'* denotes the surface area on dB{0,t) (the boundary of the ball centered at with 
radius t). For all three values of d, S{t,M.'^) = t. The initial conditions are of the form 
m(0,x) = fo{x) and ^u(0,a;) = /i(a;) for given /o, /i : M'^ ^ M. In this case, letting * denote 
convolution, 

w{t,x) = ^(5(t) * /o)(a;) + {S{t) * h){x). 

For d > 4, the fundamental solution of the wave equation is not a signed measure and so 
assumption A will not hold. 



(c) The wave equation with damping. Take L 
considered class when li < 3. Then 



+ 2a| - A on 



S{t,dy) 



^ Jo ( I a I v^2^^) 1 { I y l< t} ci?/ 

at cosh{\a\^yfi-\y\^) , 



2lT 



This also falls into the 

ff = 1, 
lid = 2, 

lid = 3. 



In these formulas, given for instance in [16] and [9], Jq and Ji are modifled Bessel functions 
of the flrst kind and of orders and 1, respectively. In these three dimensions, S{t,dy) is a 
non-negative measure. 

(d) The beam equation. In dimension d = 1, this is given by L = ^ + ^ on M. Then 
S{t,dy) = qt{y)dy, where qt{y) has Fourier transform exp(— for t > 0. The smoothness 
and integrability of qt, and hence assumption A, can be deduced from the Fourier transform 
(for example \x'^qt{x)\ < C\\d'^q/d^'^\\i). 



We now give a series representation for the solution u{t, x) of (12. ip . 
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Proposition 2.2. Let S{t,dy) be a signed measure satisfying assumption A. Suppose that V(t,y) 
and w{t,x) are bounded measurable functions on [0,T] x R"'. Define HQ{t,x) = w{t,x), and, for 
m > 0, 

Hm+i{t,x)= ds S{s,dy)V{t - s,x - y)Hm{t - s,x - y). (2.2) 
Jo Jw 

Then the integral equation h2.1\) has a unique solution satisfying sup^^^ .^.gj^d i?[|u(t, a;)p] < oo given 
by 

oo 

u{t,x) = Y,Hm{t,x) (2.3) 

m=0 

(the series converges uniformly on [0,T] x M.'^). 

Proof. We first check the convergence of the series (12. 3p . Set Mm{s) = sup^gj^d \Hm{s,z)\. Then 
Mm+i{t) < sup ds \S{s,dy)\ sup\V{r, z)\ su.p\Hmit — s,x — z)\ 

zeK'' Jo Jr'* r,z z 

< C{S,V) [ dsMUs). 
Jo 

A simple induction argument shows that Mm{s) < oo for all m,s. Gronwall's lemma (see e.g. P, 
Remark (6)]) now implies that X]m=o -^"i(^5 a;) converges uniformly on [0,T] x R*^. 

Another Gronwall argument shows the uniqueness of solutions to (12.11) . So it suffices to check 
that Ylm=o^mit,x) satisfies (12. ip . This is the case, since 



w 



„i „ oo 

{t,x)+ ds S{s,dy)V{t- s,x -y)y2Hrnit- s,x -y) 

Jo 

oo „t „ 

= Io{t,x) + y2 S{s,dy)V{t- s,x -y)H^{t- s,x -y) 

__n ^0 Jr^ 



Hmit,x). 



m=0 



Fubini's theorem, used for the first equality, applies by uniform convergence of the series and as- 
sumption A. □ 



3 Probabilistic representation 

For the probabilistic representation, we use the following additional assumption on the kernel S(t, dy) 
used in the integral equation ( 12. ip . 

Assumption B. There exists a jointly measurable process (X^, t > 0) such that for each t > 



|S((,R'')|- 

(In the case where ^(t. A) = S{t, — A), the minus sign in front of Xt is not needed). 
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Example 3.1. (a) The heat equation. In this case, one can take Xt = \fi Xq, where Xq is a 
standard X(0, 1^) random vector in . An alternative possibility is to let (XJ be a standard 
Brownian motion in W^. 

(b) The wave equation. In the three dimensional case, one can take Xj = 1 0o, where Oq is chosen 
according to the uniform probability measure on dB{0, 1). The one and two dimensional cases 
can be handled in a similar way. 

(c) The damped wave equation. Kac [T3l [H] pointed out a neat representation for solutions to 
the damped wave equation in dimension 1. Let {Na{t)) be a rate a Poisson process and define 

^ Ioi~^)^''^^^ds. If w{t,x) solves the undamped wave equation ^ — Aw = for t G M, 
X e W^, then u(t, x) = E[w{Tt, x)] solves the damped wave equation ^ + 2a^ — Au = 0, and 
with the same initial conditions. Using this, one finds that the kernel Sa{t, dy) for the damped 
equation can be written as Sait, dy) = E[S{Tt, (iy)l{rt>o} + S{—Tt, (i?/)l{^j<o}]? where S(t, dy) is 
the kernel for the wave equation. Now we satisfy assumption B by setting Xt = \Tt \ Oq, where 
00 is a uniform random variable on [—1, 1], independent of Na{t). 

(d) The beam equation. As in example (a), we can use scaling to set Xt = t^^^Xo, where Xq is 
chosen to have distribution \S{l,dy)\. 



Let X(*) = (Xf\ t > 0), i > 1, be i.i.d. copies of (X^, t > 0), and let (X(t), t > 0) be a rate 
one Poisson process independent of the (X*^*)). Let < ti < T2 < ■ ■ ■ be the jump times of {N(t)) 
and set Tq = 0. Define a process X = (X^, t > 0) as follows : 

Xt = Xo + X^'^^ forO<t <ri, 

and for i > I, 

Xt = X,^ + Xt^\ for n<t< r,+i. 

We use Px to denote a probability under which, in addition, Xq = x with probability one. 
Informally, the process X follows X*^^-* during the interval [0,ri], then follows X*^^) started at X^^ 
during [ri,r2], then X*^^) started at Xr^ during [r2,r3], etc. See Figure[T]for an illustration. 




Figure 1: A sample path of the process X and of the Poisson process {N{t)). 
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Theorem 3.2. Suppose that the kernel S{t, dy) is a non-negative measure satisfying assumptions A 
and B. Suppose w{t,x) is bounded and measurable fort < T, x G M.'^. Then {u{t,x), t <T,x & M.'^) 
defined by 



u{t,x) = e*E^ 



w 



Nit) 



t-T^^t),Xr) H \S{T,-T,^,,R'')Vit-T,,X,J 



1=1 



(3.i: 



(where, on {N{t) = 0}, the product is defined to take the value 1) is the solution of l[2.1\} . 
Proof. For t > and m > 0, let 

m 

Y{m, t) = e'l{Nit)=m}w{t - r™, XrJ JJ ["^(^^ " ^^-i' ^""Wi^ ' . 

i=l 

Then u{t,x) = X]m=o -^^[^("^' '^)]- order to show that u(t,x) is the solution of (12.11) . it suffices by 
Proposition 12.21 to show that Hm(t,x) = Ex[Y{m,t)] for all m,t,x. We prove this by induction on 
m. For m = 0, 



Ea:[Y{m,t)] = E^[eh^Nit)=o}w{t,Xo) 

= e^w{t, x)P,^{N{t) = 0} = w{t, x) = Ho{t, x). 

Now fix m > 1 and suppose by induction that Hm~i{t,x) = Erc[Y{m — for all {t,x). Set 

JFi = a{Xrj^, Ti}. Then 



E.[Y{m,t)] 
= E. 



5(ri,M'')\/(t-ri,X,Jl{.,<,}e 



n 



X E^|e* ^^l{N(t)-NiTi)=m-l}W{it - Ti) - (r„ - ri),Xr^) 
m 

H {S{{r, - n) - (t:,_i - n), M^)\/((t - n) - (r, - ri), X.J} 



X 



1=2 



i-1 



Note that, for i > 1, 

Xr^ = Xr,+^ Xi^^ + P-Tj , 

and the conditional expectation can be expressed using only the increments — ri for i > 1. Using 
the strong Markov property of {N{t)) at time ti and the independence of the families xj^^\ we deduce 
that 



E,[Y{m,t)] = E, 



S{n,R'')V{t - n, Xr,)e^n{r,<t}Y^-i{t - n, x. 

S{s,dy) 



[ ds e-'S{s,W^)e' [ 
Jo Jr" 



V{t- s,x - y)Ym-i {t- s,x -y) 



ds S{s,dy)V{t - s,x - y)H,n-i{t - s,x - y) 
Jo Vm<* 
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by the induction hypothesis and (12. 2p . This completes the proof. 



□ 



We have presented the simplest setting of the probabilistic representation, sufficient to treat our 
interest in the wave equation in dimensions d < 3 and the random potentials in the subsequent 
sections. However various extensions and variations of this representation are possible. We give a 
brief description here, leaving the details to the interested reader. 

1. For a signed kernel S(t, dy), we need to modify somewhat the representation. Write S{t, dy) = 
S^{t, dy) — S-{t, dy) for the Hahn- Jordan decomposition into a difference of non-negative mea- 
sures. Choose, if possible, subsets A(t) C M"^ so that S+(t,A(t)) = S_{t, A'^(t)) = and 



{x,t) 



[x) is measurable. (Note that this is certainly possible when S{t, dy) = qt{y)dy 



for continuous {qt{y), t > 0, y & M ).) Let Ct be a counter defined by 



i=l 



Then the argument above leads to the representation 



u{t, x) = e^Ex 



N{t) 

w{t- r^(t), (-1)^* n [ - ^-1' ■)! - ^r.) 

i=l 



where |S'(t, ■)! denotes the total variation of the measure S{t,dy). This representation then 
covers the case of the beam equation in all dimensions d > 1. 

2. If, instead of being real-valued, u(t,x) = {ui(t,x), . . . ,Un(t,x)) E M", and V(t,x) is an n x n 
matrix, so that (12. ip is in fact a system of p.d.e.'s, then the formula in Theorem 13.21 still holds, 
provided the matrix product in (13. ip is ordered according to increasing values of i. 

3. We have treated for simplicity the case of spatially homogeneous equations on W^. However, 
in principle, suitable changes should allow representations for inhomogeneous equations, or 
equations in domains with suitable boundary conditions. 

4. For any A > 0, one can replace the potential V by X~^V and use a Poisson process of rate A to 
obtain an alternative representation. For example, rewriting the heat equation (II. ip as 



du 
'di 



1 



we would get the representation 

ti(t, X 



-Au + A 



N{t) 



Wit- r^(i), X.^,^,) (-1)^* n [ - ^-1' •)! ^"'^(^ - 

i=l 



(3.2) 



where Xt starts over at the times of a rate A Poisson process. For large A, these representations, 
when using a Markovian X, become close to the classical Feynman-Kac formula. For example, 
we can further rewrite (II. ip as 



du I ^ 



u — Xu. 



(3.3) 
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Due to the term —Am in (13.31) . the Green's function e~^^pt{y) of Lm = ^ — |Am + Am gives rise 
to a factor e~^* inside the expectation in fl3.2p . which cancels the factor e^* which is outside of 
the expectation. Regarding [1 + (V^/A)] as our potential term, we find 

u{t, x) = 

Now letting A — > oo and using log(l + x) ~ x, the integrand involves a Riemann sum approxi- 
mation to the integral in the classical Feynman-Kac formula. 



w(t-r7v,(t),5,^^(,))exp 



Yl log (1 + ^ 



-1 



V{t-r^,BrJ) 



m=l 



4 Second moments for random potentials 
4.1 The random potentials 

We are now going to consider a class of linear equations driven by spatially homogeneous Gaussian 
noise F{t,x), whose covariance is formally given by 



E 



F{t,x)Fis,y) =5o{t-s)fix-y). 



In this equation, 6{-) denotes the Dirac delta function, and / : ^ M is continuous on M'^ \ {0}. 
More precisely, let V{W^~^^) be the space of Schwartz test functions (see [20j). On a given probability 
space, we define a Gaussian process F = {F{ip), ip G V{R'^~^^)) with mean zero and covariance 
functional 

E[F{ip)F{'ilj)]= [ dt [ dx f dy^{t,x)f{x-y)^{t,y). 



Since this is a covariance, it is well-known [20l Schwartz, Chap. VII, Theoreme XVII] that / must 
be symmetric and be the Fourier transform of a non-negative tempered measure /i on R'', termed 
the spectral measure : / = jF/i. In this case, F extends to a worthy martingale measure M = 
{Mt{B), t > 0, -B G i3fe(M'^)) in the sense of [2T], with covariation measure Q defined by 

Q{[0,t]x Ax B) = {M{A),MiB))t = t [ dx f dy lA{x)f{x - y)lB{y), 

and dominating measure K = Q (see [El [6]). By construction, t i-^ Mt{B) is a continuous martingale 
and 



F{^) = / Lp{t,x)M{dt,dx), 
where the stochastic integral is as defined in [21] . 



Assumption C. For each t > 0, S{t, dy) is a non-negative measure and takes values in the space of 
distributions with rapid decrease [201 Chap. VIII, §5]. Moreover, it satisfies 



rds [ fi{dO\J'S{s,-m\'<oc (4.1^ 
Jo 



and ^ 

hm/ dt [ sup \TS{rm-J'Sm)\' = 0. 

"i" Jo JR'i t<r<t+h 
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Consider the stochastic integral equation 



u{t, x) = w{t, x) + / / S{t — s, X — y)u{s,y)F{ds, dy), (4.2) 
Jo Jm.'i 

where w{t,x) is a random field satisfying appropriate conditions (see below). 

Our motivation is the case where S{t, dy) is the Green's function for a partial differential operator 
L, and the study of the stochastic p.d.e. Lu = u F, with stationary initial conditions independent of 
F. This s. p.d.e. can be recast into this integral form with w{t,x) being the solution of Lw = with 
the same initial conditions as u{t,x). In this context, {w{t,-), Mt{-)) is stationary in x, or, more 
precisely, has property (S) of Dalang [6, Definition 5.1]. 

The stochastic integral in (14. 2 p needs defining. If S{s, y) is a smooth function, as in the heat 
equation, then we can use the stochastic integral with respect to a worthy martingale measure 
introduced in [21]. In this case, (14.21) has a unique solution provided w{t,x) is a predictable process 
such that sup«y ^^-^^d E[w'^{t,x)] < oo. 

If S{s, ■) is a singular measure, as in the case of the 3-dimensional wave equation that we are 
particularly interested in, then we use the integral introduced in Dalang [6]. We briefiy describe 
his construction, that uses an approximation to the identity. Choose ip G C^(R'^) with ip > 0, the 
support of ip is contained in the unit ball of M"^ and f^^ ip{x)dx = 1. For £ > 1, set ^ipe^x) = C'^tp^lx), 
so that — i> 5o as £ — > oo. The stochastic integral in (14. 2 p is the L^-limit of the usual stochastic 
integrals 



/ / Se{t - s,x - y)u{s,y)F{ds,dy). 

Jo jRd 



where Si{t, x) is the convolution / S{t, dy)ipe{x — y). While studying the s. p.d.e. Lu = uF as above, 
this convergence was established in [6], and the same arguments show existence and uniqueness of a 
solution to (14. 2 p provided w{t, x) has the property (S) of [HI Definition 5.1] and supj<2. E[w'^{t, 0)] < oo 
for all T > 0. Details for this can be found in [Bj Section 5]. Assumption C (in particular, the fact that 
S{t, dy) is non-negative) is also used in the definition of the stochastic integral. In the terminology 
of [9j, u(t,x) is a random field solution of (14. 2p . that is defined for every t and x (as opposed to a 
function-valued solution, defined only for all t and almost all x, that would not be adequate for our 
purposes). 

In fact, it is shown in [6] that (14. ip is even a necessary condition for (14. 2 p to have a solution 
satisfying supj<y ^gjgd -^[^^(t, x)] < oo. 

In the cases of the heat and wave equations, [6] gives equivalent conditions to (14. ip involving only 
H or the function / in the covariance structure. 



4.2 The series representation 

In this subsection, we work under assumptions A and C. We assume that w{t,x) has the properties 
indicated in subsection 14 . 1 1 that ensure that the stochastic integral in (14.20 is well-defined and ensure 
existence and uniqueness of a random field solution to this integral equation. 

We shall show that there is a series representation for the solution u{t,x) of (14. 2p . analogous 
to (12. 3p . but with the deterministic integral replaced by the stochastic integral as in (14. 2p . Define 
Io(t,x) = w(t,x), and, for m>0, 

Im+i{t,x)= j / S{t- s,x -y)Im{s,y)F{ds,dy). (4.3) 
Jo Jr'^ 



11 



Proposition 4.1. Suppose that w{t,x) is bounded and measurable for t < T, x & W^. Then the 
series 



u{t,x) = ^Im{t,x) (4.4) 

m=0 

converges in uniformly over {t,x) G [0,T] x and is the unique solution to Ili4.2\ ). 
Proof. We first check the L^-convergence of the series in fl4.4p . Set 

M„(t) = snp E[Im{t,xy]. 



By P, Theorem 2], 

MUt)< [ dsM„,.i{s) [ fi{dO\J'S{t~s,-m\\ 

Jo jRd 

By (14. ip and [SI [71 Lemma 15], we conclude that 

oo 

MmisY^^ <oo, 

m=0 

which establishes the L^-convergence of the series. Set Un(t,x) = J2m=o^rnit,x). Then Un(t,x) 
u{t,x) in L^, and by [SI Theorem 2], as n — oo, 

S(t - s,x - y)un{s,y)F{ds,dy) ^ / / S{t - s,x - y)u{x,y)F{ds,dy). 

Jo ^Rd 

Therefore, 

w{t,x)+ / / S{t — s, X — y)u{s,y)F{ds, dy) 



lim / / S{t - s,x - y)un{s,y)F{ds,dy) 

V Jo Ju'i 

lim i Io{t,x) + / S{t- s,x -y)Im{x,y)F{ds,dy) 



n+1 



= lim Imit,x) 

m=0 

= u{t,x), 

showing that u{t,x) solves fl4.2p . □ 

The successive terms in (14. 4p are orthogonal in L^, that is E[Im{t,x)Im'{s,y)] = whenever 
m 7^ m'. The series is therefore a chaos expansion for the noise F. The orthogonality can be checked 
by induction on m and m', using the fact that the covariance between Im and Im' reduces, as in (14. 5 p 
below, to an expression involving the covariance between Im-i and Im'-i- 
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4.3 The probabilistic representation of second moments 

In this subsection, we work under assumptions A, B, and C. We make the same assumptions on 
w(t,x) as in subsection 14.21 

Let {N{t), t > 0) be a rate one Poisson process. Using two independent i.i.d. famihes {X^^'^\ i > 
1) and (X*-*'^'', i > 1), construct, as in Section[3l two processes = (X/, t > 0) and = (X^^, t > 
0) which renew themselves at the same set of jump times Tj of the process N, and which start, under 
Pxi,x2y ^1 ^2 respectively. See Figure [2] for an illustration. 




H 

To = 



H ^ 



n r2 



T-3 



Figure 2: A sample path of the processes X^ and X^ 



Theorem 4.2. Let u{t,x) be the solution of given in Proposition 4jJ_ Then 
E[u{t, x)u{t, y)] = e^E^^y w (t - rjv(t) , X^^^^^ ) w (^t - r^(t), X^^^^J 

Nit) 

1=1 

Proof. Observe that by Proposition 14.11 

oo oo oo 

E[u{t,x)u{t,y)] = J2Y1 E{Im{t,x)Im'{t,y)] = E[Im{t,x)Im{t,y 



m=Q m'=0 



m=0 



using the orthogonality of the terms in the series. For m > 1, using the smoothed kernels Se{t,x) 
defined earlier, we have 

E[Irr,{t,x)Im{t,y)] 

= E S{t ~ s,x — z)Im-i{s, z)F{ds,dz) / / S{t — s^y ~ z)Im-i{s, z)F{ds,dz) 

Jo JRi Jo Jr'^ 

= Urn E / / Si{t - s,x - z)Lm-i{s,z)F{ds,dz) / SfXt - s,y - z)Im~i{s, z)F{ds,dz) 

L^O JR-i Jo JRd 

= lim ds dzi dz2 Se{t - s,x - Zi)Si>{t - s,y - Z2)f{zi - Z2)E[I^_i{s, Zi)lm-i{s, Z2)] 
Jo Jr^ Jr^ 

= ds S{t-s,x-dzi) S{t - s,y - dz2)f{zi - Z2)E[Im-i{s,zi)Im-i{s,Z2)], (4.5) 
Jo Jr'' Vr<* 



13 



where we have used the Lebesgue Differentiation Theorem [221 Chapter 7, Exercise 2] in the final 
step. We shall now show by induction that 



E[Imit,x)Imit,y)] = J{m,t,x,y), m > 0, 



(4.6) 



where 



J{m,t,x,y) = e^E, 



x,y 



kN{t)=m}W {t - Tm, X^^) W {t - T^, X^J) 

m 

X{{S{n-n-,X?f{xl^-xl)} 



X 



i=l 



For m = 0, 



J(0,t,x,y) = e'w{t,x)w{t,y)P,^y{N{t) = 0} = E[Io{t, x)Ioit, 



We suppose now that (14. 6 p holds for m — 1. By the Markov property at ri, arguing as in Theorem 
13.21 we have, choosing J^i = (j{ti, Xl_^, X^_^}, 



J{m,t,x,y) 



E.,y[l{r,<t}e^'f{Xl^-X^JS{T,,R''f 

X e^^'''^^E^,y[l{N{t)-N{T^)=m-l}W{t - Tm,XUw{t - T^, X^J 
m 

X Ul^{r.-r,^i,^rfiXl^-Xl)\T, 

i=2 

E.,y [lir,<t}e^'fiX,, - Xl)SinXfJim -l,t- t^.X\.X\ 



ds / S{s,x — dzi) I S{s,y — dz2)f{zi — Z2)J{m — l,t — s, zi, Z2). 
Jr'' 



The conclusion now follows from (14.51) and the induction hypothesis. 



□ 



Remark 4.3. By multiplying the integral formulas (14.21) for u{t,x) and u{t,y) and taking expecta- 
tions, one expects formally the integral equation 



E[u{t,x)u{t,y)] 
= w(t, x)w{t, y) + 



j f f S{t - s,x - dzi)S{t - s,y - dz2)f{zi - Z2)E[u{s,Zi)u{s,Z2)]. 
Jo Jr'^ Jr"^ 



This new integral equation on M^o! is 

of the same form as (12. ip . This leads to an alternative derivation, 
by applying Theorem 13.21 of the representation for second moments given above. However, we have 
used the argument above as it will generalize to higher moments. 



5 Moments of order n 

In this subsection, we work again under assumptions A, B, an C. In addition to the assumptions 
on w{t,x) made in subsection 14. 2^ we assume that sup^^j. E[\wit, x)\'p] < 00, for all T,p> 0, which 
ensures that the solutions have finite p-th moments. 
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In the case where u{t, x) solves a first order equation driven by the Gaussian noise F, written in the 
form du/dt = Lu + uF, then a formal calculation suggests that the n-th moment m(t, Xi, . . . , x„) = 
E[u{t, Xi) . . . u(t, Xn)] should satisfy 



^ = L^^,...,xr,rn + ^ f{xi - Xj) 



(this formula is proved for discrete space in [3l Section II. 3]). Here, Ivxi,...,x„ stands for the sum of 
the operator L applied to each variable x,. The equation for m is again of the same potential type 
considered in section [21 and can be recast as an integral equation using a multiple product kernel 
constructed out of the kernel S(t, dy) for L. Theorem 13.21 then leads to a probabilistic representation 
for m. This argument does not seem to apply for second order equations or directly for integral 
equations. However, as we shall now explain, it is possible to find a representation, analogous to the 
one for second moments, that holds for the higher moments of the general integral equation (14. 2p . 

We start with an informal discussion of the representation for higher moments. The second 
moments were given in terms of a pair of processes, both of which were renewed at the times of 
a single Poisson process A^. The situation for the n-th moment is somewhat analogous. Instead of 
two processes, we will use n processes X^, . . . ,X"'. For each pair of indices p = {pi,p2}, we create 
a Poisson process Xt(p). The renewal times of the the process X* will be the union of the Poisson 
times arising from the processes Nt{p), such that the index i is contained in the pair of indices p. 

More precisely we let Vn denote the set of unordered pairs from = {1, . . . , n} and for p G P„, 
we write p = {pi,p2}, with pi < p2. Note that card (Vn) = n{n — l)/2. Let {N.{p), p G Vn) 
be independent rate one Poisson processes. For A ^ Vn let Nt{A) = X]pGA^i(p)- This defines a 
Poisson random measure such that for fixed A, {Nt{A), t > 0) is a Poisson process with intensity 
card(A). Let ai < (72 < ■ ■ ■ be the jump times of {N^iVn), t > 0), and i?* = {R\,R2} be the pair 
corresponding to time (Xj. Two possible representations of this Poisson random measure are shown 
in Figure [31 

For i G Cn, let P*-^^ C Vn be the set of pairs that contain i, so that card(P*^^^) = n — 1. Let 
rf < r| < ■ • ■ be the jump times of (Xi(pW), t > 0). We write Nt{i) instead of Xt(pW). Note that 

5^ Ntip) = N,{Vn) = IY1 ^*(^)- 

We now define the motion process needed. For £ G £„ and i > 0, let (Xf^'\ t > 0) be i.i.d. copies 
of the process (X^) defined before Example 13. 1[ Set 



Xo^ + Xf«, 0<t<rf, 



y\^e -t- ^j_^£ , ' i ^ ^ 



This motion process is illustrated in Figure [H 

It will be useful to define X/ for certain t < 0. For given (ti,xi), . . . , (t„, under the measure 
P{h,x,),...,{tr.,x„) we set 

X/ = X^i°) for -U<t< 0. 

Finally we set Tq = —ti. The following theorem gives a formula for the n-th moments, and it is the 
main result of this section. 
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{3,4}.. 
{2,4}.. 

{2,3} 
{1,4} 
{1,3} 
{1,2} 



H 1 1 1 h- 

(Ji (72 CTS Cr4 CTS 



H 1 1 1 h- 

(Ji 0-2 0-4 0-5 



Figure 3: Two equivalent representations of the Poisson random measure {Nt{-)): the top rep- 
resentation is simply the superposition of the Poisson processes {Nt{p)), p G P„; in the bottom 
representation, two elements of £„ are joined at time CTj if they constitute the pair i?*. 



Theorem 5.1. The n-th moments are given by 
E[u{t,xi) ■ ■ ■u{t,Xn)] 

= gtn(n-l)/2 ^ 



(5.1) 



NtiVn) 

n /(^?-^?) 

i=l 



^(0,xi),..,(0,x„) 

Nt{e) 



e 



The proof of this theorem requires some preliminaries. Let x), m > 0, be as defined in fl4.3l) . 
For < s < t, set 



Im+l{s, t, x) 



S{t -r,x- y)Im{r,y)F{dr,dy), 



JR'i 



SO that Imit,t,x) = Imit,x) foY 771 > 1. Foi 171 = and < s < t, we set Io{s,t,x) = Io(t,t,x) 
w{t, x). Let 



I{s; {rrii, ti,Xi),i = 1, . . . ,n) = E 



^rrti ("5, ^j, Xi 



,i=l 



for mj > 0, s < min(ti, . . . , t„), Xi G M°', i = 1, . . . , n. We begin by giving an inductive expression for 
this expectation. 

Lemma 5.2. Suppose mi + ■ ■ ■ + m„ = m. 
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H h 

tI 



To 



To 



Ti To 



Figure 4: Illustration of the motion processes in the case where n = 4 and Xq 



X 



1,...,4. 



(a) If m = 0, then 



I{s; (0, ti,Xi),i = l,...,n) = Y[ w{ti 



(h) If m > 1, then 

I{s; {mi,ti,Xi),i = l,...,n) 



E 



X E 



JR'i 



i=l 



dr / ^(tp, - r, rfyi) / S{tp^ - r, dy2)f{xp^ - Vi - Xp., + 7/2) 



fce£n\p 



i=l 



(5.2) 



Proof. Part (a) follows immediately from the definitions. For part (b), if m = 1, then n — 1 of the rrii 
are equal to and so n — 1 of the ImX^^ti^Xi) are deterministic. The one Imi{s,ti,Xi) with irti = 1 
is a martingale with mean zero, implying that J(s; {■mi,ti,Xi),i = 1,. . . ,n) = 0. The expression in 
formula (15.21) is also equal to since there is no p e Vn such that rrip-^ ■ rrip^ > 0. 

If m > 2, we distinguish two cases. The first case is where all but one of the rrii are zero. In this 
case, /(s; {mi,ti,Xi),i = 1, . . . ,n) and expression (15.21) vanish, for the same reasons as in the case 
m = 1. We now consider the second case, in which there is at least one p G P„ with rrip-^ ■ rrip^ > 0. 

Using the smoothed kernels Sg = ipi* S, as in section 14.11 we define 



Si>{t -r,x- y)Imir, y)F{dr, dy). 



For fixed (tj,Xj), s ^ I^_{s,ti,Xi) is a martingale, and according to [211 Thm. 2.5], if mj > and 
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rrij > 0, then the mutual variation process of ll^^{-,ti,Xi) and ll^,{-,tj, Xj) is 



dr / dyi Si{ti -r,Xi- yi) / dy2 Se{tj - r, xj - 1/2) 
X f{yi - ?/2)/m,-i(r, r, ?/i)/„^_i(r, r, 7/2)- 

We now apply Ito's formula [H Theorem 5.10] to the function /(ai, . . . , a„) = ai ■ ■ ■ a„ and the n 
martingales ll^.{-,ti, Xi), i = 1, . . . ,n. Note that 



9V 
da! 



and 



^ fce{l,...,n}\{i,i} 



JJ ttfc if i 7^ j. 



The stochastic integrals terms given by Ito's formula have mean zero, because the Jmi(-,-,-) have 
bounded moments of all orders, so taking expectations we reach 



E 



Y[ll^{s,ti,Xi 



i=l 



dr / dyi Se{tp^ - r, Xp, - yi) dy2 Se{tp^ - r, Xp., - 2/2) 



X f{yi - y2)E 



Imp^-i{r,r,yi)Imp^_i{r,r,y2) JJ I^^^{r,tk,Xk) 

k€C„\p 



(5.3) 



The variables Im^ tk, Xk) and /^^ (r, tk, Xk) are both bounded in for all p and continuous in in 
(r, Xfc), so that they are continuous in in the variables (r, x^). This implies that the expectation in 
(15.31) is continuous in (r, xi, . . . , x„). Using the change of variables zi = Xp^ —yi and Z2 = Xp^ —y2, we 
let £ 00 in (15. 3p . The left-hand side converges to I{s; {rrii, ti,Xi),i = 1, . . . , n) and the right-hand 
side converges to formula (15. 2p . completing the proof. □ 

Define 



J{t; {mi,ti,Xi),i = 1, . . . ,n) 



gin(n-l)/2 



(ti,Xi),...,(t„,Xri) 



i(miH hm„) 



l{Nt{e)=me 



i=l 



X I I I I ^{T- - T, 

ieCn i=l 



i-l! ■ 



The next aim is to show that these expectations satisfy a similar inductive formula. 
Lemma 5.3. Suppose mi + ■ ■ ■ + m„ = m. 

(a) If m = 0, then 

n 

J{t; (0, ti,Xi),i = 1, . . . ,n) = J]^ w(t + ti, Xi). 
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(b) If m>l, then J(t; {mi,ti,Xi),i G £„) is equal to 

y2 ds S{tp^ + s, dyi) / S{tp^ + s, dy2)f{xp^ - Vi - + 1/2) 

X J{t- s; {rup^ - 1,0, Xp^ - yi), i = 1,2; {rrii, s + te, Xe), ^ e \p). (5.4) 



Proof. Part (a) follows immediately from the definitions. For part (b), in the case that only one of 
the rrii are non-zero then J(t; {m(,,ti,Xi),i G £„) = since P{ti,xi),...,{t„,x„){Nt{i) = rrii, £ e = 0, 
and formula (15. 4p is also equal to since there is no p G P„ with m^^ ■ rrip^ > 0. 

We now suppose that m > 2 and that there is at least one p E Vn with rrip^ ■ rrip^ > 0. In this 
case {Nt{i) = mi, i G £„} C {ai < t}, and we are going to use the Markov property of NtiVn) at 

time (Ti. Indeed, choosing jFi = o{ai,R ,Xa^ ,Xa^}, we may rewrite J{t; {mp,te, xe), i G £„) as 



-E'(ti,a;i),...,(t„,x„) 



m 

X -E'(ti,xi),...,(t„,x„)[l{Aft(^)-Ar<,,(£)=m^-, <?G£„\p}n{Aft(i;)-Af<,,(<?)=mf-l, tep} H /(-^(^!' " ^(^^^) 



X n ns( 

£6£„\p i=l 



I _ I 



nns( 



o-l' 



1,2, while for 



Note that at time cxi, on {R^ = p}, the processes X^' start afresh from X^ 

i E Cn\ p, X^ has seen no jump from —ti to ai, that is for di + t^ units of time. Using the strong 
Markov property at cxi, the conditional expectation above multiplied by e(*~'^i)"("~^)/^ is equal to 

J(t - (Ti; (m^^ - 1,0, X^;), i = 1,2; (m^, di £ G Cn\p)- 

Therefore, J(t; {mi,ti,xi),i G £„) is equal to 

pG7^,i:mpj-mp2>0 ^Sp 

X J(t - cTi; (mp, - 1,0,X^;), i = 1, 2; (m^, cJi + t^, x^), £g£„\p)]. 

The variable ai is exponential with mean 2/ (n(n— 1)) and the variable R is independent and uniformly 
distributed over Taking the expectation over cxi, i?\X^; we reach ([El. □ 



Proof of Theorem \5.1\ We note that it suffices to prove, when ti > t, that 

lit; (mi,ti,xi), . . . , (m„,t„,x„)) = J(t; (mi,ti - . . . , (m„,)f:„ - (5.5) 
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oo oo 



Indeed, in this case, by Proposition 14.11 

E[u{t,Xi) ■ ■ -uit^Xn)] = ^ ■ ■ ■ E{Im,{t,t,Xi) ■ ■ ■ Ira„it,t,Xn)) 

mi=0 m„=0 
oo oo 

= X] ■ ■ ■ 5Z . . . , (m„, t, Xn)) 

mi=0 mn=0 
oo oo 

= 5Z -^^l^^; ("^i>0,xi), . . . , (m„,0,x„)), 

mi=0 m,„=0 

and this is equal to the expression in (15.11) . 

Let m = mi + - ■ ■+mn- We are going to prove (15. 5p by induction on m. If m = 0, then (15.51) follows 
from Lemma [5^ a) and Lemma IS^T a). since both sides of ( 15. 5p are equal to w{ti,xi) ■ ■ -witn^Xn)- 
Now assume inductively that (15.51) holds for m — 1 > 0. By Lemma [5.3( b). 

J{t] {rrii, ti -t,xe), i e Cn) 

= ds S{tp^ -t + s,dyi) S{tp,, ~t + s,dy2) 

X /(a^pi -yi- + 2/2) (5.6) 
X J{t- s; {nip^ - 1,0, - yi), i = 1,2; (m^, s + - t.x^), le Cn- p)- 

By the induction hypothesis, the last factor J{t — s; . . .) is equal to 

I{t - s; {mp^ - 1, t - s, Xp^ - yi), i = 1, 2; (m^, t^, xt), i e Cn\p)- 

Now use the change of variables r = t — s and Lemma 15.2( b) to see that the right-hand side of ( 15. 6p 
is equal to I(t; (m^,t^,x^), i G £„). This completes the proof. □ 



Remark 5.4. The intuition behind equality (15.51) is the following. Suppose n = 4 and consider space- 
time positions (ti, xi), . . . , {t^, X4), as in Figure [5l The quantity I(t; (m^, tg, xg), £ = 1, ... 4) is the 
expected product of iterated integrals, where the left-most integral is up to time t < min(ti, . . . ,^4) 
and the order of the iterated integrals are mi, . . . , 777,4. 

On the other hand, time s for the Poisson random measure runs in the opposite direction as in 
the s.p.d.e. (see Figure |5]). In the quantity J(t, {rrii, ti — t,xe), £ = 1, ... 4), the process starts at 
negative time t — tg, and there are no Poisson pairs during negative time. During the time interval 
s = to s = t, the number of Poisson pairs containing xe is set to 777^. With these constraints, 
I(t; {mi, ti, xi), £ = 1, ... 4) = J(t, {nil, ti — t,xi), £ = 1, ... 4) as stated in (15. 5p . 
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